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SUMMARY 


In this article, we consider the application of three popular domain decomposition methods to 
Lagrange-type nonconforming finite element discretizations of scalar, self-adjoint, second order 
elliptic equations. The additive Schwarz method of Dryja and Widlund, the vertex space method of 
Smith, and the balancing method of Mandel applied to nonconforming elements are shown to 
converge at a rate no worse than their applications to the standard conforming piecewise linear 
Galerkin discretization. Essentially, the theory for the nonconforming elements is inherited from the 
existing theory for the conforming elements with only modest modification by constructing an 
isomorphism between the nonconforming finite element space and a space of continuous piecewise 
linear functions. 


INTRODUCTION 

We consider the convergence properties of domain decomposition methods applied to 
Lagrange-type nonconforming finite element discretizations ol scalar, self-adjoint, second order 
elliptic problems. An isomorphism between the nonconforming finite element space with the 
natural norm induced by the elliptic problem and a conforming piecewise linear space with the 
i/ 1 - semi norm is constructed. Using the isomorphism, we are able to apply the existing analysis of 
domain decomposition methods for conforming elements to nonconforming elements with only 
modest modifications. As examples of this technique, we show that the operators arising in three 
popular domain decomposition methods, specifically, the additive Schwarz method of Dryja and 
Widlund [1], the vertex space method of Smith [2], and balancing method ol Mandel [3], applied to 
nonconforming finite elements have condition numbers that satisfy the same bounds as the ones 
given in [4] and [5] for conforming finite elements. 

The same technique was used in [6] and [7] to analyze the rate ol convergence of balancing 
domain decomposition and the standard additive Schwarz method for the dual-variable mixed finite 
element formulation. Moreover, as a corollary of the analysis of Smith’s method for the 
nonconforming spaces presented in this paper, we have a new bound for Smith s method applied to 
mixed finite elements. 

♦This work was supported in part by the National Science Foundation under Grant No. DMS-9112847. 
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After the research for this paper was completed, the author was made aware of some related 
work done concurrently by Sarkis [8]. In particular, the isomorphism used herein was independently 
suggested by Sarkis for linear nonconforming elements. In [8], Sarkis constructs and analyzes 
special coarse spaces such that when the overlapping additive Schwarz method is applied, the 
condition number of the resulting operator is bounded by a constant times (1 + \og(H/h))(l + H/S ) 
in both two and three dimensions. Here H and h are the characteristic sizes of the subdomains and 
mesh, respectively, and S is a measure of the overlap of subdomains. The notable characteristic of 
Sarkis’ bound is that the constant is independent of jumps in the coefficients across subdomain 
boundaries. If the techniques of this paper were used to derive bounds that were independent of 
the jumps in coefficients, the resulting bound would include one log factor in two dimensions using 
[1, 9], but two logs in three dimensions using [5, 10, 11]. 

The remainder of this paper is divided into six sections. In the next section, we set some 
notation, formulate the nonconforming problem, and construct an equivalent representation in 
terms of the nodal values. In Section 3, we construct an isomorphism between the nonconforming 
space and a continuous space of piecewise linear functions. The isomorphism is used in Section 4 to 
analyze the rate of convergence of the Dryja-Widlund additive Schwarz method. In the last three 
sections, we consider the substructuring methods of Smith and Mandel applied to the 
nonconforming problem. 


PRELIMINARIES 


We consider the following self-adjoint, uniformly elliptic problem for p on the polygonal domain 
C IR n , n = 2,3, with boundary dfl: 


—V AVp = f in fi, p = 0 on <9fl, 


(i) 


where A is a uniformly positive definite, bounded, symmetric second order tensor, and / G L 2 (f2). 
The uniform ellipticity of (1) implies the existence of positive constants c„,c* such that the 
following bound holds: 


< £M(*x < cTt ve € nr.Vx'G n. 


( 2 ) 


In order to set a length scale, we assume that the diameter of is one. We introduce a two level 
quasi-regular triangulation of 0: a division first into subdomains {H;} *L X with diameter 0(H ) , and 
a refinement of the first into elements with diameter 0(h). Following [12], define the scaled Sobolev 
norms 

1 II .lO ii i • *"> ..T I....O 


1 1 , 0 , 


IN 


H 2 




u 


1/2, an, “ Mi/ 2 , an, + ^ IMIo 


i.an, 


where 


IMIo, n, = / |«(x)| 2 </x, IMIo, an, = / |u(s)| 2 rf«, 

J n, Jan, 


t — s \ n 


(It ds. 


Let .V(fl) be a finite dimensional nonconforming finite element space of Lagrange-type defined 
subordinate to the triangulation T that vanishes at all degrees of freedom on <9fb Since Af( fi) is of 
Lagrange-type, the elements in Af(Vl) may be expressed in terms of a nodal basis, and we may 
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identify an element in Af(il) with the values it attains at the nodal points. For convenience, we 
assume that the subdomains and the elements are triangular in two dimensions or tetrahedral in 
three dimensions. Extensions to other shape regular decompositions are straightforward. 

We consider the problem of finding p/, 6 Af(Cl) such that 

<l{Ph,<]h) = f fqudx W//, € (3) 

Jn 

where d is the generalized Dirichlet form: 

d{ph,<lh) = dn(Pk,qu), da>(ph,qh)= X AVp h ■ ¥q h dx. 

t£T, TCfi' JT 

We now introduce several conventions used in this paper. In this paper, we shall only be 
concerned with the solution of this finite dimensional problem, and will henceforth drop the “h” 
subscript. 

Having defined a parent finite element space of functions A'(fl) with a nodal basis and a set 
fP C fl, we will simply write A*(f V) for the restriction of ,T(fi) to iT, i.e. 

A'(SY) = {0|o* | <f> 6 A'(fi)}. 

By an abuse of notation, we consider an element. <j> G A'(fi') also to be an element of A'(fl) by 
setting <j) to zero at all nodes outside of Q f . 

We will write Q\ ~ Q 2 if two quadratic forms Q 1 and Q 2 with the same domain T> are 
equivalent , i.e. if there exists constants ci,c 2 > 0 such that 

ciQi(<t>, 4>) < Q2{4 >i < f>) < c 2 Qi{<f>, <j > ), V<£ e T>. 


In what follows, C will be used to denote a generic constant that may not be the same from one 
line to the next. This constant, as well as the constants involved in the equivalence of quadratic 
forms, will always be independent of h and //, but can depend 011 the constants in (2), the shape 
regularity of the subdomains, the degree of the nonconforming finite elements, and the regularity of 
the triangulation. 

To conclude this section, we prove a lemma that provides an equivalent quadratic form for d(*, •) 
in terms of the nodal degrees of freedom. The proof of this lemma was suggested by Joseph Pasciak 
in the context of the mixed finite methods considered in [6, 7]. 

Lemma 1 Let f}' C Q be the union of elements of T . And let A(x) — a(;r)A(.x) ; where a is a 
positive } piecewise constant function with value a T on r E T . Then for every p £ .A/"(f2') ; 

< M/MO - X X (/'("•) - (4) 

r ET, nodes : 

r C Q f n, , 7ij £ T 


The constants that appear in the definition of the equivalence do not depend on the constants in (2), 
but rather on constants that arise when A is replaced by A. 


Proof. The local kernel of d T (-, ■) in W(r) is exactly the constant functions on r since for 

p € V(r) 

d T (p,q) = 0 W/ G Af(r), <=> Vp = 0. 
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Figure 1; Refinement of the ‘2D P-1 element and a partial refinement ol the 3D P-1 element. 

Hence, (c? T (-, •)) 1 ^ 2 is a norm on A^(t)/IR. Since all norms are equivalent on finite dimensional 
spaces, we see that 

d T (p,p) ~ a T |r| 1_2/n (P( n d ~ P( r h)) 2 1 

nodes : 

Tli , Tlj € T 

by a simple scaling argument. The proof is completed by summing over the elements of T in ft'. 


□ 


A CONFORMING EQUIVALENCE 


In this section, we construct a conforming space that is isomorphic to V(ft) using the techniques 
in [6, 7] and recall some basic properties about the isomorphism. 

Given an element t £ T, let Tj. be a subtriangulation ol t such that the vertices of the 
subtriangulation include the vertices of t and the nodal points in t pertaining to the degrees of 
freedom of V(t). Every element in the new triangulation should have at least one vertex that 
corresponds to a nodal point ol V(t). Moreover, the subtriangulations should be constructed in 
such a way that the union of subtriangulations gives rise to a refined quasi-regular triangulation of 
ft which we denote by 

f = U f r. 

reT 

A vertex of T will be called primary if it was a nodal point corresponding to a degree of freedom of 
V(ft); otherwise, we call the vertex secondary. We say that two vertices of the triangulation T are 
adjacent if there exists an edge of f connecting the vertices. An example of the subtriangulation of 
the P-1 element that has nodal degrees of freedom at the center of its edges (faces) is given in 
Figure 1. 
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Let (4(f)) denote the space of continuous piecewise linear functions subordinate to the. 
triangulation T that vanish on 50. For 0' C 0, a union of elements, define £4( O') by restriction, i.e. 

£4(0') = {u\n> | « 6 £4(0)} • 

Since the functions in £4(0') are naturally parameterized by the values they attain at the 
vertices, we can define a pseudo-interpolation operator 2°' into £4(0') for any function <j> defined at 
the primary vertices contained in 0' by 


f 0, if x e 50' n 50; 


<f>(x), if x is a primary vertex not in 50' D 50; 


p'm = 


The average of all adjacent primary vertices on the boundary 
of O', if x is a secondary vertex in 50' \ 50; 

The average of all adjacent primary vertices, if x is a secondary 
vertex in the interior of O'; 


(5) 


The continuous piecewise linear interpolant of the above vertex 
values, if x is not a vertex of T. 

Since 2°' is well defined for any function defined at. the primary vertices, by an abuse of notation, 
we can understand 2°' both as a map from ./V’(O') into £4(0') and a map from £4(0') into itself. 

For any 0' that is the union of elements in T, let £2/ t ( O') C £4(0') denote the range of T ; that is, 

££/.(0') = {^-2 n V/, f /€A£(0)}. 

We now prove that 2°' : -> £4(0') preservesj.lie norm induced by the bilinear form •) 

on W(O') and the H 1 - seminorm on £4(0'). Since f 1 ' is a bijection between A£(0') and £4(0') by 
construction, this proves that A^O') and £4(0') are isomorphic. 

Theorem 2 Let 0' C 0 be the union of dements. Then for all p G A£(0'), 

Pli.n'- 

Proof. This proof is an expanded version of the proof given in [7]. Recall that for <j) € £4(0'), 

| E M 1 ’ 27 " E ( ? ) 

r £ 'f vertices : 

rCtl' W.w, €r 


By virtue of Lemma 1 and Equation (7), it is enough to show that 


E M 1_27n E (pM-pK)) 2 - 

t £ T, nodes : 

tCD' n t ,Tij e t 


E_ £ 

T £ T vertices : 

7 -CO' Vi,Vj€T 


((2 fi 'p)(^) - (2 ft 'p)(^)) 2 - 


Since vertices of T t contain the nodal points ol r and p — 2° p at these points, we have 


E (p( n ») -pM) 2 ^ C E 

nodes : T67V 

Tli ,Tlj 6 r 


E ((?'/-) w - (2”'p) to))’ . 

vertices : 

, Vj £ T 


(8) 
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where the constant is controlled by the regularity of the subtriangulation. Hence, by summing over 
the elements of T in H', we conclude that the right hand side of (8) dominates the left hand side. 

To prove that the left hand side dominates the right, we note that the differences in the right 
hand side are of three types: the difference at two primary vertices, the difference at two secondary 
vertices, and the difference at a primary and a secondary vertex. Since p and 2 .p agree at primary 
vertices of T, the difference at two primary vertices occurs as a term in the left hand side. For two 
secondary vertices v\, in an element r £ T containing a primary vertex v p , we see that 

((P'p) („,) - (Pp) (, 2 )) 2 < 2 ((2"p) („,) - (2“V) M) 2 + 2 ((i°'p) M - (F P ) („,))’ . 

Hence, it is enough to bound the difference at a secondary and primary vertex by terms in the left 
hand side of (8). 

Let v n+I be a secondary vertex with adjacent primary vertices v 1} . . . , v n7 and let p 3 = p(v 3 ). 
Noting that for j — 1, . . . , n 

P ° P) (c) = Pj » (^V) K+l ) = - {? r p) K ) = 

U J = 1 11 J=1 

we see that 


((^’p) {Vn+l) ~ {fp) (V ,)) 2 


1 


n* 


it,(pj - Pi)] 


n 



Ti 




by the Cauchy-Schwarz inequality. The proof is completed by summing over all triangles of T . The 
number of such terms, and hence the constant in the bound, is controlled since the regularity of the 
mesh implies that there is an a priori maximum number of adjacent elements that can share a 
secondary point. □ 

Using the techniques in the proof of Theorem 2, the following lemma is easy to prove. 


Lemma 3 There exists a constant C depending only on the regularity of the triangulation T and 
the degree of the nonconforming space such that for any IT C the union of elements of T, 

I^Vkn- < C\<f>\ kt m \/<p € £/,.(«'), k = 0,1. (9) 


THE DRYJA-WIDLUND ADDITIVE SCHWARZ METHOD 

The presentation in this section and the next follows the treatment of Schwarz methods given by 
Dryja and Widlund in [4j. We concentrate only on the additive Schwarz methods with exact solves. 
The convergence rate of the multiplicative Schwarz method may be estimated in terms of the same 
quantities (see [13]) and is easily worked out. Extensions to inexact solves are likewise direct. 

Recall that the additive Schwarz method with exact solves for (3) is completely determined by a 
decomposition of the finite element space V(D) = Ao + Mi + . . . + Mm- For each subspace M t , 
define an operator Pi : — > M, by 

d(PiP, <l) = d(p, q) V<y G Mi. (10) 
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The additive Schwarz algorithm with exact solves for (3) involves the solution of 

M M 

pp = f \ p = E / = £/<, (ii) 

i— U i=0 

where /, G Afi is defined by 

= [ fqdx Wq € M- 
Jn 

Abstract bounds on the condition number of P have been derived in terms of two quantities, Co 
and the spectral radius of £ , which we now define. Let Co be a constant such that for every p G Af 
there exists a representation p = Yiio Pi with p x G Afi satisfying 

M 

- C 0 d(p t p). (12) 

i=0 

Let p(£) denote the spectral radius of £ = {t^}, the matrix of strengthened Cauchy-Schwarz 
constants; that is, is the smallest constant for which 

\d(Pi,Pj ) | < €ijd(pi t pi)id(pj,pj)* Vpi € Mi, Vpj € Mj, i,j > 1. (13) 

The next theorem, due to Dryja and Widlund [14], bounds the condition number of the additive 
Schwarz method in terms of Co and />(£): 

Theorem 4 The eigenvalues and the condition number n(P) of P satisfy 

Amin (P) > C£\ X max (P) < (p(£) + 1), K(P) < Co(p(£) + 1). (14) 

To construct the decomposition of A^(fi) to be used in our application of the additive Schwarz 
algorithm for nonconforming elements, we first create an overlapping decomposition of the domain 
ft by extending each subdomain ft, to a larger region ft- which is also the union of elements of T . 
We characterize the extent of the overlap of the partition {ft-}^ x by 6, where 

5 = min dist(<9ft t - \ 5ft, 5ft' \ 5ft). 

The decomposition {ft-}-^ gives rise to a natural decomposition of Af(£l) by letting Af t C Af(Sl) 
denote the set of functions that vanish at all nodes in the closure of (ft \ ft'). In order to provide a 
mechanism for global exchange of information between subdomains so as to enhance the rate of 
convergence, we also use a low dimensional space defined by 

Mo = {p € M( n) | p = <£ e U H ( fl)}, 

where is nodal interpolation into M{£l), and f///(fl) is the space of continuous functions that are 
linear on each subdomain ft,. Note that the subspac.es tor the nonconforming space are exactly the 
nodal interpolants of the standard decomposition of the conforming space f//,(ft), namely, 
t4(ft) fl i/o(ft'). 

In the following lemma we recall the crux of the proof due to Dryja and Widlund (Theorem 3 of 
[4]) that the Schwarz method applied to the conforming Galerkin discretization has a condition 
number that is 0(1 -f ( H/8 )). 
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Lemma 5 For every € U h (Cl), there exists a decomposition <f = <f>* with 4> o € U w (ft), 

<j>i £ £4(ft) H 1 < i < M and a constant C independent of h, H , and 8, such that 

- C ( } + y) 

1=0 

We now show that the application of the Schwarz method to the nonconforming space converges 
at the same rate. 

Theorem 6 The condition number k(P) of the additive Schwarz operator P defined by (11) induced 
by the decomposition V(ft) = Afo + ■ • • + Am of the nonconforming finite element space satisfies 

4P) < c (l + j) ■ 

The constant C is independent of h, 8, and II . 

Proof. The verification that the largest eigenvalue of P is bounded by a constant is standard. 
Since d{p„ Pj ) = 0 for p, € A f it Pj € Mj with fi' n ft'- = 0, P may be written as the sum of an a 
priori bounded number of disjoint projections. Since projections have unit norm, a constant bound 
on the largest eigenvalue of P is immediate. See, e.g., Lemma 3.1 of [2]. 

For P € A /'(ft), let (7 Q P )i denote the decomposition of 2°p € U h (Sl) arising in Lemma 5, and set 
Pi = {{iPpfi). It is easy to check that p, € A f and p - E^oP>- Using Theorem 2 and Lemma 3, 
we see that for * = 0, , M, 

d(p„ P .) < cii n ((i n ri.)i;.n < ci(2 n l .).i;,„. 

Summing and applying Lemma 5 and Theorem 2, we conclude that 

MM / ff\ ~ / H\ 

Edte.pt) <C£ KPphlla <C[ 1 + j) iz°pl?.n < C (l + j) d(p,p). 

i= 0 i = 0 

Hence, C 0 in (12) is bounded by C (1 + H/8). An application of Theorem 4 completes the proof. □ 

SUBSTRUCTURING DOMAIN DECOMPOSITION 


The remaining two methods considered in this paper are domain decomposition methods applied 
to a reduced problem involving only the degrees of freedom on the internal interfaces of subdomains 
p _ 3ft; \ 5ft. Following [4], we recall the construction of the reduced problem. Since AZ'(ft) is 
of Lagrange-type, we may associate with functions p, q € V(ft) the vectors of values they attain at 
the nodes. Let x and y denote the vectors of nodal values of p and q, respectively, and x l ’,y ' the 
subvectors of degrees of freedom in ft,. Let D (,) denote the local stiffness matrix arising from 
dn.(-,-), and let D denote the global stiffness matrix, i.e. 


x <'> T D<V> x T Dv= Z D'%" = <!(,,, q). 

t =l M 
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For each subdomain, we can partition the degrees of freedom into two sets, the ones related to 
nodes on the boundary of ft; denoted Xg\ and the ones corresponding to nodes in the interior of ft, 
denoted x^\ Such a partitioning induces a partitioning of D ^ given by 


x M T DWyW = 


The interior unknowns of each subdomain may be eliminated in terms of the boundary unknowns. 
The resulting matrix, 5", is the Schur complement with respect to the interface unknowns defined by 



xlSy B = Y, : 

Af 


B 


s® 


Vb 


where 


= d§ B - Dfj(D';n 


(0\-i n (0 

U IB- 


It will be convenient to work with the bilinear forms induced by S and S^'\ and so we define 

s(p,(j) = xlSy B , *(/»,?) = x fl S (,) j4 ) . 

For a function p G Af(fi), we note that unlike conforming spaces, the restriction of p to the 
interfaces, pjn, is not solely determined by the nodal values on T since A^fi) is nonconforming. 
Hence, we are careful to understand Af(T) as a subset of Af(fi) parameterized by the nodal values 
on T consisting of the discrete harmonic extension of the nodal values to the interior of the 
subdomains. Specifically, if p G AT(T) has the vector of nodal values x ^ on d Q t , then p|n t is the 
function associated with the vector of nodal values {x\ l \x$) T where x\'^ satisfies 


u ii x i 


_n(0 T (d 
U IB X B * 


A linear functional g is easily constructed such that finding p G AZ^T) satisfying 


s(p. i) = sM V, € V(0 


(16) 


is essentially equivalent to (3). 

We now construct a conforming space of functions that is isomorphic to A r (F) with the norm 
induced by the bilinear form s(-,-). Let U/ L ( T) denote the restriction of Uh(Q) to U^l 1 dQ t . Since 
functions in U/ l ( T) vanish on dfi (because functions in do), functions in Uh(T) can be 

parameterized in the natural nodal basis by the values they attain at the vertices of T in I\ 
Analogous to (5), for T f the union of edges (and faces in 3D) in the triangulation T and <f> a 
function defined at the primary vertices in P, define a pseudo-interpolant T <f> G Uh{T f ) by 

f 0, if x G P fl dQ] 


<f>(x ), if x is a primary vertex not in V fl dCl] 


T<t>{x) 


The average of all adjacent primary vertices on T f if x is a 
secondary vertex on P; 


(17) 


The continuous piecewise linear interpolant of the above vertex 
values, if x is not a vertex of T. 
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Note that if T = dW, then iF'f = (2^)|an< for all <j> in jV’(ft') that agree with <{> at the nodal 
degrees of freedom of dfl'. 

Since jT is well defined for any function defined at primary vertices, by an abuse of notation, we 
can understand if both as a map from A^L') into Un (F) and a map from Un{ H) into Uh{T'). We 
denote the range of jF by 

tMH = «?V)ir# € t/,,(r')). 

The equivalences in the following lemma are a combination of the standard trace theorem and an 
extension theorem for U h (dQ t ). In particular, the proof of this lemma given in [6] shows that the 
space Uk(Fli) is rich enough to inherit the Extension Theorem of Widlund [15] from £//,(fl,)- 

Lemma 7 For <f> € Uh{dfli), 

II ^11 1/2,30, ^ inf IMkn.-, kli/2,3n, ^ inf l^li.n,- (18) 

<t> e u h (Q t ) 4> e u u {Q,j 

4>\on t - <t> ^|sfi, - $ 


Additionally, there exists a constant C independent of mesh, parameters such that 

|r> n $ M o, < c% m W e Unidn,), k = o, 1 / 2 . (19) 

The following theorem plays the role of Theorem 2 for the interface problem. 

Theorem 8 For all p € A/^T), 


»i(p,p) ~ |2 on 7 , li/2 1 ?n,‘ 


(20) 


Proof. By a direct computation followed by an application of Theorem 2 and Lemma 7 noting 
that Uk(Fli) = ^(A/^fi,-)), we see that 


s t 


i(?,P ) 


inf d Qi (p, p) 
P6 AT(n.) 

Pjan, = Pt 


~ _ inf l^'pipn, - l^ n 'Pli/2,3fi- 
pjzAT(n.) 

Pjcin, = V 


□ 


SMITH’S VERTEX SPACE METHOD 


Smith’s vertex space method [2] is an additive Schwarz method applied to the interface problem 
(16). The decomposition of Af( P) is constructed slightly differently in two and three dimensions. In 
both cases, we first partition T into overlapping subsets based on its decomposition as the boundary 
of subdomains. In two dimensions, for each vertex Vj of T, let denote the set of points on T 
that are less than a distance S from V 3 . For each edge Ei of T, let Tf' denote the interior of the 
edge E{. In three dimensions, for each vertex V 3 , each edge E{, and each face Fk of T, define T S J as 
above, let Tf* denote the interior of the face F k , and let Tf' denote the set of all points in strips of 
width 8 on all faces which share the common edge Ei. 

Understanding the set of faces to be empty in two dimensions, the decomposition of T into 
subsets induces a decomposition of V(T) by considering 

Mr) = £ ), 

GelH,Ei,Vj,F k } 
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where for G G {E{ y Vj, Af(Tf) C A/^r) are those functions that vanish at all nodal points on T 
that are outside of the set Tf , and Af(T ^ ) C Af(T) are those functions that are the nodal 
interpolant of the restriction to T of continuous functions that are linear on each subdomain and 
vanish on <90. 

The following lemma is the crux of the analysis of Smith’s method by Dryja and Widlund [4] for 
conforming elements. 

Lemma 9 For every <f> € Uf t { T), there exists a decomposition 

<t> = Y fa 

G€{H,E„V Jt F k } 

with cf>„ G <fc € U h (Tf) = Un(T) n H' 0 (Tf) for G G {£„ V„ F k ) such that 

M M 

E El*ol?/ W o, SCIl+kffl/JD’EKtan,. (21) 

Ge{H y E ly V 3 ,F k } t=l i = 1 

The constant C is independent of the choice of (j>, and the mesh parameters h } H } and 6. 

Let Pp : A/^T) — ► Af(T) denote the additive Schwarz operator defined by (10) with the bilinear 
form d(*, ■) replaced by the interface form $(•, •) and the decomposition of Af(fl) replaced by the 
decomposition of Af(T) described above. We now prove that the condition number of A/r for the 
nonconforming space has the same bound given in [4] for the similar operator for the conforming 
finite element space. 

Theorem 10 The condition number of the additive Schumrz operator P p for Smith’s decomposition 
for the nonconforming finite element discretization satisfies 

K(Pr) < C((l + log (H/h)) 2 . (22) 


The constant C is independent of the mesh parameters h } H , and 6. 


Proof As in the proof of Theorem G, Pr may be written as the sum of an a priori bounded 
number of disjoint projections, and so the largest eigenvalue of Pp is bounded by a constant. 

To bound the smallest eigenvalue, we also proceed as in the proof of Theorem 6. For p G Af(T), 
set p G = Jp((^p) G ), G G [H, Ei, V u F k }, where if is interpolation at the nodes on T into A/*(r) 
and (x p)c is the decomposition of Tp G UffiT) that arises in Lemma 9. Since Tp and p agree at 
the nodal degrees of freedom of A/"(r), and 

X(T") = 1^\U h (V)), U(T r ;) = l^(V, l (T r n) VG e 

it is easy to check that 

p= Y p°- 

Ge{H,E„V J ,F k } 

Working one subdomain at a time and using Theorem 8 and Lemma 3, we see that for G = H and 
for G G [E{, Vj, Fk} such that Tf fl dQ, ^ 0 we have 


s,(K,K) < ClFVl?^ = G|f’ (1 '((F, l ) 0 )i; /2iMi < CKfrtol?^. 


(23) 
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Assume that we can prove that there exists a constant independent of h, H and 8 such that 

M M 

£ PpfiAon, < C ■£ ^■rWn.m. V,. e -V(r). ( 24 ) 

1 = 1 * = 1 

Then by summing (23) over subdomains and subspaces, noting that Si(pa,pa) = 0 if fl = 0, 
and applying Lemma 9, Equation (24), and Theorem 8, we see that 

M M 

^2 s(pg,pg ) = 5Z £ s i(PG,PG) < C (1 + log {H/8)) 2 ^2 l^pli/ 2 ,an< 

Ge{H,E„Vj,F k ) Ge{H,E„Vj,F k }i=l i=1 

M 

< C (1 + \og(H/8)) 2 J2\l da 'P\ 2 i/2,dn, < C(1 + log (H/8)) 2 s(p,p). 

i = 1 

The proof of the condition number bound now follows from an application of Theorem 4, and we 
are only left to verify (24). ■ ' 

Define a pseudo-interpolant, : Af{Cl) —* f//, ( 0 ) by (5), noting that the boundary of D \ T is 
dfl U r. Using the techniques in the proof of Theorem 2, it is easy to show that there exists a 
constant C\ depending only on the regularity of the mesh and the degree of the nonconforming 
space such thaU 

M M 

£ I^W.n. < C. £ Vj> € W(«). 

t=l 1=1 

By Lemma 7, for each p € A/"(T) there exists an extension p E 6 A/"(f!) that agrees with p at the 
nodal points on T such that 

PVlU < C\T™'p\\ ri>mi i = 1 , . . . , M. • 

Combining these results alter another application of Lemma 7 with <j> = dfip, we conclude that 
M M M 

£i?pi;/wo, <c£if ! vi;,< 1 . <c£i2"'p E i;,n. 

i=l 1 = 1 1 = 1 

which verifies (24). □ 

In [6], the interface form arising from the discretization by mixed finite elements of (1) was 
shown to satisfy Theorem 8 with A'(L) replaced by the appropriate space of interelement 
multipliers. Hence, the proof given above is applicable to discretization by mixed finite elements, 
and we arrive at the following corollary. 

Corollary 11 The application of Smith’s decomposition method to the dual-variable mixed finite 
element formulation discussed in [6] results in an operator whose condition number grows at worst 
like <?((! + log {H/8)) 2 ). 


BALANCING DOMAIN DECOMPOSITION 


As the final domain decomposition method considered in this paper, we investigate the balancing 
domain decomposition method of Mandel [3] applied to nonconforming finite elements. The method 
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involves the iterative solution (usually by conjugate gradients) of (16) preconditioned by the 
balancing preconditioner described in Algorithm 1 below. Each iteration involves the solution of a 
local problem with Dirichlet data, a local problem with Neumann data, and a “coarse-grid” 
problem to propagate information globally and to insure the consistency of the Neumann problem. 
The theory and practical performance of balancing domain decomposition for the standard 
conforming Galerkin finite element method and mixed finite element method are the subjects of [5] 
and [6], respectively. As in previous sections, we will deduce the convergence theory for the 
nonconforming spaces from the conforming theory in [5] using the isomorphism introduced in the 
fifth section of this paper. 

One remarkable property of balancing domain decomposition is that the bound on the condition 
number of the preconditioned operator is independent of jumps in coefficients across subdomains. 
Specifically, let the tensor A in (1) be written as A{x) = a(;r)A(x), where a is a positive function 
that is piecewise constant with constant value a, on ft;. The uniform ellipticity then implies that 
there exists positive constants c*,c* such that 

c.a,t T Z < ?A{x)t < caii T ( € IR n ,Vx € f C (25) 

The bound on the condition number of the operator that arises in balancing domain decomposition 
will depend on c* and c* but will be independent of a, and c* and c* in (2). 

Following Mandel’s original exposition in [3], we now recall the balancing preconditioner in terms 
of matrices. A equivalent variational presentation is given in [6]. By an abuse of notation, we use 
the same symbol to denote an element in Af(T) and its associated vector of values attained at the 
nodal degrees of freedom. 

The balancing preconditioner is parameterized by two sets of matrices, a set of weighting 
matrices {Wi}^ and a set of kernel generators {Zi)^L x . The weighting matrices 
Wi : Af(dfli) —> Af(dQi) are chosen such that they form a decomposition of unity on A^F), i.e. 

M 

Y,N t w % N?p = p VpeW(n, 

1=1 

where N{ denotes the canonical inclusion mapping TV; : A/"(3ft;) — > M(T) by extending elements of 
A/^dft;) by zero at all other degrees of freedom. A prescription for the weighting matrices that 
guarantees a convergence bound independent of coefficient jumps between subdomains is given in 
Lemma 12 below. For each subdomain ft,, let n; — dim(Af(3ft;)), and select an n, x m x matrix Z x of 
full column rank with 0 < m, < n,, such that 

KerS, C Range Z;, i = 1, . . . , M. (26) 

For the scalar, second order, elliptic problems we consider in this paper, KerS*; is empty if there is 
Dirichlet data imposed on any part of 5ft, Pi 3ft, otherwise it is the set of functions that have the 
same value at all the nodes on 3ft,. From the kernel generators, we construct a “coarse space”, 

Mh C jV^r), defined by 

M 

M H = { P e Af(T) : p = £ W V,z,ze Rang eZ,j. 

1 = 1 

We say that q € Af(T) is balanced if it is orthogonal to A///; that is, 

Zj W?Nfq = 0, i = 1 , . . . , M. (27) 
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The process of replacing r by a balanced q = r — Sw, w £ A///, will be called balancing and involves 
solving a “coarse grid problem” over the space A T/ described in (28) and (29) below. 

The action of the balancing preconditioner is defined by the following algorithm. 

Algorithm 1 Given r £ AT T ) , compute Mfff^r as follows. Balance the original residual by solving 


the auxiliary problem for unknown vectors \ 3 £ IR mj , 

M 

zj W? Nj {r- SJ2 NjWjZjkj) = 0, i = 1, . . . , M (28) 

j = l 

and set 

M 

q = r — S ^2 NjWjZjXj, qi = W?Njq, i = (29) 

j= 1 

Find any solution G ^(<9(7,) for each of the local problems 

SiUi = q iy i = 1, . . , A/, (30) 

balance the residual by solving the auxiliary problem for p 3 G JR m > } 

M 

%W?HT(r - «£ + ZjW)) = 0. i-1 M. (31) 

i=i 

and set 

M 

^ = p^u,U/4 (32) 

i= 1 


If some rrij = 0, then Z 3 as well as the block unknowns \i 3 and A^ are void and the j- th block 
equation is taken out of (28) and (31). 

In [3], it was proven that Algorithm 1 implements a well defined operator that is symmetric and 
positive definite. An abstract bound on the condition number of S was also given. We will use 
the following lemma proven in ~[5] to determine a bound on the condition number of the 
preconditioned system for the application to the nonconforming discretization. 

Lemma 12 For subdomain f7 l; define the weighting map W l as multiplication of the nodal values 
such that for any pi G 9^(317*) and each nodal point v G <9Q t; 

(Wipi)(u) = — =5 P«-(0- (33) 

2_, Q i 

{j:n t £dQ 3 } 


Assume that there exists a number R so that 


-s,(NjN.p„ 


Nj N t p t ) < — R S{(pi,pi) 

on 


(34) 


for all i, j = 1, . . . , M and all p x G Af(dQ{) that are orthogonal to the range of Z{. Then there exists 
a constant C not dependent on h } H or R } so that the condition number of the 

preconditioned system satisfies 

< C R. 
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Since most nonconforming methods, like the mixed finite elements considered in [6], do not have 
nodal degrees of freedom at vertices in two dimensions and vertices and edges in three dimensions, 
the analysis of [6] is directly applicable to the nonconforming case. However, in keeping with the 
philosophy of this paper, we will allow vertex and edge degrees of freedom and will deduce the 
general nonconforming theory from the conforming theory. The following definitions and lemma 
from [5] will provide the essential conforming theory that we need for the nonconforming case. 

Definition 1 Any vertex, edge, and , in the 3D case, face, of the interfaces between subdomains 
{H;} will be called a glob. A glob is understood to be relatively open; for example, an edge does not 
contain its endpoints. We will also identify a glob with the set of the degrees of freedom associated 
with it. The set of all globs will be denoted by Q . 

Definition 2 For a glob G, define the selection operator Eg * U\ t (T) — ► Uh(T) as follows: for 
<f> G Uk(T), Eq 4> is the unique function in U} t (T) that has the same values as <f> on the degrees of 
freedom in G , and all other degrees of freedom of Eo<f> are zero. 

Note that the union of all globs disjoint ly cover the set of all degrees of freedom of [/^(T), and 
the mappings Eg are projections that induce a decomposition of unity on Uk(T), E e g = ’■ 

GeG 

Lemma 13 There exists a constant C such that for all p* G Af(dDi) that are orthogonal to the 
range of Z x and for all globs G G D dflj 

£ C(1 + 

Proof The proposition follows from Lemmas 3.7, 4.6 and 5.1 of [5]. □ 

We now prove a bound on the condition number for the preconditioned system in balancing 
domain decomposition for nonconforming elements. 

Theorem 14 The interface operator S preconditioned by the balancing preconditioner Mb a i defined 
in Algorithm 1 with weighting maps defined in (33) has a condition number k(M^S) satisfying the 
bound 

K(M^S)<C(l + \o g (H/h))> (35) 

in both two and three dimensions with the constant C independent of h, H, and a t . 

Proof By tracing back the dependence on the coefficients, it is easy to prove the following 
refinement of Theorem 8: 

s,{p,I>) =si o,|r 9n 7 '|i /2iani Vp € (36) 

with equivalence constants that no longer depending on c * and c* in (2), but only on c* and c* in 
(25) - 

For pi 6 we may decompose Nf NiPi as the nodal interpolant of globs: 
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Let |p,||. = 5,(p,,pi). Considering those p, that are orthogonal to the range of Z,, working one glob 
at a time, and using (36), (19), Lemma 13 and (36) in that order, we have 

\I^f {E a ^'Pi)\% < a jC\T >Qj EcT )Q '[>i\i/2,dn ] 5: a j^'\EG^ 0Q 'Pi\] /2,dQ, (37) 

< a i C(l + lo6(///A)) 2 |5 0 % i |? /2ian< 

< — C(1 + log (H/h)) 2 \pi\ 2 s .. 

OCi 

By the construction of the decomposition, there is an a priori maximum number of globs that 
intersect dH, H dClj. Summing over such globs, we conclude that 

Sj{Nf Nifi, Nj NiPi) < —C(l + log (H/h)) 2 Si {pi,pi). 

J J a, 

The proof is completed by appealing to the bound in Lemma 12. □ 
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